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1 . 1 PURPOSE 

This  report  seeks  to  describe  experiments  performed  in  September 
1976  with  an  existing  LWIR  sensor  at  the  AMOS  Observatory.  The  goal  of 
these  measurements  is  an  estimate  of  wavelength  scaling  of  angle- of- arrival 
fluctuations  for  radiant  energy  propagating  via  slant  paths  through  the  whole 
atmosphere.  The  overall  object  of  this  study  is  to  provide  information  rela- 
tive to  defining  physical  limits  imposed  by  turbulent  atmosphere  upon  imaging 
and  collimating  systems,  principally  at  LWIR  wavelengths. 

1.2  BACKGROUND 


The  angle  of  arrival  experiment  is  part  of  a wider  program  of  mea- 
surements begun  under  ARPA  sponsorship,  through  their  agents,  RADC,  in 
July  1975.  In  its  entirety,  the  IR  Atmospheric  Measurements  Program 
(Contract  No.  F30602- 75- C- 0235)  also  included;  (2)  real-time  LWIR  at- 
mospheric absorptance  determination;  (3)  a demonstration  of  the  effective- 
ness of  long  period  signal  integration  for  distinguishing  deep- space  satellite 
signatures  from  gradients  associated  with  intervening  radiant  atmosphere; 

(4)  measurement  of  wavelength  scaling  of  atmosphere  affected  point  image 
spread  function;  and  (5)  measurement  of  the  spectral  broadening  of  'mono- 
chromatic' (ultra  stable  CO2  laser)  radiation  produced  by  propagation  through 
the  entire  turbulent  atmosphere.  The  equipments  for  accomplishing  items 
(2)  and  (3)  (which  were  designed  and  fabricated  at  AERL  as  part  of  this  pro- 
gram) were  left  in  place  at  AMOS  and  have  been  used  subsequently  in  support 
of  satellite  observations.  A brief  description  of  this  hardware  and  some 
preliminary  results  are  presented  in  Appendix  A and  References  (1)  and  (2). 

In  May  1976  reevaluation  of  then  current  needs  and  risks  produced 
a decision  to  cancel  items  (4)  and  (5).  The  remaining  experiments  were  re- 
scheduled for  earlier  completion. 

1.  Avco  Everett  Research  Laboratory  Interim  Report;  IR  Atmospheric 
Measurements,  Contract  No.  F30602-75-C-0235  (December  1976). 

2.  Young,  L.A.,  IR  Radiation  in  a Real  Tropical  Atmosphere,  p.77  (May  1977). 


1.3  NEED 


The  projected  performance  of  surveillance  and  other  systems  which 
depend  upon  LWIR  radiative  transfer  over  slant  paths  through  the  entire 

atmosphere  are  presently  based  upon  theoretical  extension  of  visual  tur- 

(3) 

bulence  data  obtained  from  observations  of  stars  and  upon  monochromatic 
measurements  along  horizontal  paths.  The  detailed  configuration  of  large 
adaptive- optic  collectors  and  projectors  for  imaging  through  atmosphere  will 
hinge  upon  the  scale  of  refractive  structure  and  no  truly  representative  data 
for  confirming  the  utility  of  thecretmal  scaling  is  available.  Indeed  there  is 
some  observational  evidence  which  conflicts  with  prediction.  (See  Appendix 
B.) 

1.4  MEASURING  CONCEPT 

The  measuring  principle  for  angle- of- arrival  represents  an  extension 

of  the  familiar  Hartmann  test  for  determining  the  imaging  performance  of 

(4) 

large  telescopes  in  the  field.  In  1904  Hartmann  described  a procedure, 
using  bright  stars,  which  is  less  subjective  than  the  older  knife-edge  (Foucault 
testing)  observation.  He  suggested  masking  the  entrance  of  the  telescope  to 
be  characterized  with  an  opaque  screen  pierced  with  a carefully  mapped  pat- 
tern of  smallish  holes.  When  light  from  a distant  point  source  is  imaged  by 
the  masked  optic  a photographic  plate,  positioned  some  distance  from  focus, 
records  a fascimile  of  the  hole  pattern  as  distorted  by  imperfections  of  the 
optic  and  random  atmospheric  refraction  structure.  Any  measurable  differ- 
ences between  the  resulting  photograph  and  the  screen  pattern  are  then  used 
to  trace  rays  from  each  hole  and  to  characterize,  thereby,  the  performance 
of  the  optic.  Atmospheric  effects  are  a nuisance  which  limit  the  precision 
of  these  tests.  To  average  out  refractive  disturbances,  produced  by  wind 
driven  eddies,  long  exposures  (typically  greater  than  ten  seconds)  are  made 
for  optical  testing. 

3.  Tatar skij,  V.I.,  Effects  of  Turbulent  Atmosphere  on  Wave  Propagation 
NtlS  YT68- 50464  (1971). 

4.  Hartmann,  J.  , Zeit.  f.  Instrum.  24  (1904). 
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We,  on  the  other  hand,  sought  to  characterize  the  turbulent  atmo- 
sphere. The  quality  of  the  optic,  as  long  as  it  remained  constant,  was  of 
little  importance.  This  was  accomplished  by  making,  in  effect,  a great 
many  short  exposure  records  and  using  the  difference  between  records  to 
determine  atmosphere  induced  fluctuations  in  the  apparent  direction  of 
arrival  of  star  light  across  the  aperture. 

1.5  HARDWARE  DESCRIPTION 

Instead  of  cine  film  we  used  an  array  of  detectors  in  an  existing 
radiometer  (AMTA)  and  a correspondingly  coarse  pattern  of  holes  as  shown 
in  Figure  1. 

Sixteen  of  the  twenty-five  AMTA  detectors  are  used  for  this  experi- 
ment. Each  spot  is  aligned  so  that  it  average  position  corresponds  to  the 
common  intersection  of  four  detectors;  by  comparing  the  signal  histories 
for  each  member  of  a quadrature  set  radial  components,  amplitude  and 
direction,  of  spot  motion  may  be  derived.  The  correlation  between  signal 
histories  for  detectors  monitoring  different  spots  is  a measure  of  the  degree 
to  which  spot  motions  are  tethered.  If  the  atmospheric  refractive  (or  ray 
tilting)  structure  is  large  compared  to  the  telescope  entrance  size,  the 
spots  ought  to  move  as  a tethered  ensemble.  If  signals  between  even  close 
spaced  spots  are  uncorrelated,  refractance  structure  is  small  compared  to 
fifteen  inches.  Repeating  the  observation  for  different  colors  yields  data 
for  calculating  structure  function  scaling  with  wavelength. 

Table  1 lists  the  wavelength  options  exercised  during  this  series  of 
observations.  Atmospheric  transmission  and  source  temperature  effects 
are  included.  The  three  principal  objects  observed  were  Mercury,  Venus 
and  a • Bootis  with  effective  temperatures  of  560°K,  240°K  and  4900°K, 
respectively. 

Because  our  measurements  extend  in  various  discrete  bands  to  the 
LWIR  (1.  6 to  21.  5 jum),  thermal  emission  from  an  ambient  temperature 
mask  in  the  entrance  aperture  would  have  overwhelmed  radiation  from  even 
the  most  intense  astronomical  sources.  Therefore,  a small  screen  (detailed 
in  Figure  2)  refrigerated  to  40°K,  was  positioned  within  the  AMTA  radiometer 
vacuum  dewar  as  shown  in  Figure  3 near  an  image  of  the  telescope  pupil. 
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TABLE  1. 


SPECTRAL  OPTIONS 


Source 

Temp 

240°K 

560°K 

4900 

°K 

Filter  # 

^ eff 

AX 

^ eff 

AX 

^eff 

AX 

1 

2 

3.  9jLim 

4.8 

0. 42fj  m 
0.41 

3. 78jjm 

4.  8 

0. 67pm 
0.42 

3. 66pm 

0.  81pm 

3 

8.74 

0.91 

8.70 

0.98 

4 

11.63 

2.62 

11.47 

2.44 

5 

11. 04 

3.55 

10.  51 

3.81 

6 

18.74 

2.98 

18.63 

2.84 

9 

4.  26 

1.02 

3.75 

1.10 

1.6 

0.62 

Light  imaged  by  the  telescope  is  brought  to  focus  again,  after  being  spatially 
quantized  by  the  screen,  approximately  midway  between  the  screen  and  the 
detector  array;  see  Figure  13.  Focus  was  relocated  by  moving  the  radiom- 
eter reimaging  mirror  as  shown  in  Figure  2.  It  was  also  necessary  to  move 
the  screen  and  the  pupil  image  so  the  former  would  mask  interfering  self- 
radiance of  the  warm  scan  mirror  as  shown  in  Figure  3 and  this  was  accom- 
plished by  increasing  the  power  of  the  KBr  vacuum  window.  All  these 
modifications  were  designed  to  be  quickly  reversible  so  the  sensor  could  be 
restored,  without  degraded  performance,  to  its  normal  measuring  configu- 
ration. The  holes  in  the  screen  were  sized  to  project  approximately  eight 
inch  diameter  subapertures  into  the  telescope  entrance.  Figure  4 illustrates 
the  results  of  computations  relating  the  diffracted  hole  projections  in  the 
array  plane  to  effective  entrance  diameter  for  representative  wavelengths. 
Detector  spacing  and  the  longest  wavelength  constrained  the  hole  size.  For 
each  wavelength  shown  there  is  an  effective  entrance  diameter  below  and 
above  which  the  illuminated  spot  on  the  detectors  increases  in  size  due  to 
physical  and  geometric  optical  effects,  respectively.  [Note:  There  are  two 
curves  for  visible  light  representing  absolutely  still  and  normally  turbulent 
(1.  5 arcsecond  "seeing")  conditions.  The  Hartmann  screens  for  telescope 
evaluation  have  two  inch  holes  and  test  photographic  plates,  which  are  also 
exposed  approximately  six  inches  from  focus,  typically  have  . 035  to  . 040  inch 
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Figure  2 Optical  Layout 
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Figure  3 AMTA  Dewar 


diameter  spots  when  seeing  conditions  are  'good.  ' Good  seeing  means  no 
observable  twinkle  (unaided  vision)  above  45°  elevation.  ] The  ordinal 
values  in  Figure  4 are  radius  for  80%  encircled  energy. 

Signals  from  all  twenty-five  AMTA  detectors  are  independently  and 
simultaneously  recorded  on  analog  magnetic  tape  with  AC  bandwidth  extend- 
ing from  0.  2 Hz  to  4.4  kHz  for  each.  Steady  state  values  are  not  recorded 
and  signal  level  is  adjusted  so  that  detector  noise  (sky  radiance  produced) 
may  be  retrieved  from  tape.  Certain  housekeeping  information  like  Universal 
Time,  spectral  option  (filter)  number  and  telescope  gimbal  error  are  also 
recorded.  Telescope  line  of  sight  wander  during  the  experiment  would,  of 
course,  produce  correlated  motion  of  the  subaperture  projections.  To  per- 
mit this  effect  to  be  unwound  from  the  data,  provisions  for  recording  gimbal 
shaft-angle  encoder  error  signals  were  made. 

The  AMOS  telescopes  are  mounted  in  high  speed  (2°  per  second) 
three-axis  (azimuth,  polar  and  declination  rotations)  gimbals.  Shaft  angle 
values  for  any  two  active  components  are  read  with  23  bit  (0.  15  arcseconds) 
resolution  1,  000  times  per  second  and  compared  with  continuously  recomputed 
desired  values;  the  differences  are  converted  to  analog  voltages.  These, 
called  'mount  error  signals,  ' were  patched  to  the  analog  sensor  recorder 
for  this  experiment. 

1.6  MEASUREMENT  CONSTRAINTS 

Two  identifiable  sources  of  difficulty  were  anticipated:  (a)  mount 
servo  instability,  and  (b)  sensor  refrigerator  "thumping”  which  produces 
microphonic  signals  correlated  between  all  detectors.  The  first  named 
produces  oscillatory  or  random  (wind  buffeting  related)  fine  motion  about 
the  desired  line  of  sight.  If  the  telescope  steering  servo  gains  are  set 
'properly'  the  mount  will  oscillate  between  five  and  ten  Hertz  with  an  ampli- 
tude ranging  between  twenty  and  'a  few  tenths'  arcseconds.  If  the  gains  are 
backed  off  to  decisively  preclude  oscillation  the  mount  becomes  compliant 
enough  to  respond  to  wind  disturbance  torques.  Instability  is  generally  con- 
fined to  the  declination  axis;  resonant  or  phase  shift  producing  elements  in 
the  feedback  loop  had  not  been  unambiguously  identified  in  September.  The 
effect  is  generally  adjusted  and  left  so  as  to  produce  little  or  no  degradation 
for  photometric  or  radiometric  observations.  We  hoped  to  circumvent  the 
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problem  by  using  the  error  signal  record  and/or  selecting  favorable  dome 
slot  orientations  and  windless  nights  for  observations. 

Sensor  microphonics,  i.  e.  , response  to  refrigerator  component 
acceleration,  affects  high  frequency  (80  and  219  Hz)  interference  which  we 
assumed  would  be  in  phase  for  all  detectors.  The  effect  is  most  pronounced 
for  detectors  far  off  the  optical  axis  and  is  principally  due  to  small  vibration 
indiced  rocking  of  the  scan  mirror  and  to  movement  of  the  detector  array 
itself.  To  reduce  microphonic  effects  the  scan  mirror  was  blocked.  This 
did  reduce  the  effect  but  at  insupportable  cost;  we  needed  to  steer  the  beam 
with  the  scan  mirror  to  establish  initial  alignment.  Instead  the  problem  was 
reduced  during  data  processing  by  notching  out  at  the  principle  offending 
frequencies. 

Figure  5a/b  relates  tilt  measuring  sensitivity  to  source  brightness. 
These  plots  are  based  upon  calculated  spot  motion  and  upon  noise  equivalent 
flux  measurements  made  two  years  ago  for  filters  #5  and  #6.  Taking  1/40 
wave  of  tilt  with  20%  measuring  uncertainty  as  a goal  leaves  only  a dozen 
candidate  source  objects.  Objects  had  (a)  to  be  brighter  than  m^  = -6; 

(b)  to  be  unresolvable  by  the  eight  inch  subapertures  at  the  wavelengths  of 
interest;  (c)  be  far  enough  from  the  sun  to  avoid  heating  the  telescope's 
mirrors  (>  12°)  and  finally;  (d)  be  > 2°  from  the  moon  to  assure  the  possi- 
bility of  visual  guidance.  We  found  that  sensor  NEFD  performance  had 
degraded  by  a factor  of  four  to  five  from  predictions  based  upon  the  old 
measurements  and,  in  consequence,  restricted  our  candidate  list  still 
further  eliminating  all  but  the  terrestrial  planets  for  LWIR  observations 
and  Q • Boo  for  near  IR.  Figures  6 and  7 give  geometric  observational 
particulars  for  the  planets  and  maximum  allowable  image  size.  Mars  is 
sufficiently  bright  only  for  long  wave  measurements.  Mercury  is  bright 
enough  at  all  wavelengths  during  our  September  1976  measurement  slot  but 
comes  too  close  to  the  sun  between  the  14th  and  24th  days.  Venus  is  bright 
enough  also,  but  is  of  marginal  size  for  the  short  wave  observations. 


5.  AERL:  Mission  Instruction  and  Operation  Plan  #018  (attachment  to 
A/OP-309)  (July  1976). 
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Figure  5 Detectable  Wavefront  Tilt 
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Figure  7 Diffraction  Blur  Size 


The  observations  had  to  be  terminated  whether  or  not  data  was 
gathered  by  16  September  to  permit  restoring  the  sensor  in  preparation 
for  satellite  radiometric  missions  scheduled  shortly  thereafter. 
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2.  0 CALIBRATION 


■f  1 


To  relate  recorded  observational  data  to  angle  of  arrival,  a 'point' 
radiance  source,  fitted  with  micronneter  cross-slides,  was  positioned  in 
the  telescope  Cassegrain  image  surface  to  establish  the  relation  between 
detector  signals  and  image  displacement.  Tilt  in  microradians  is  given 
by  = 6s/F;  where  F and  6s  are  effective  focal  length  and  off-axis  image 
displacement,  respectively.  The  cross-slide  had  sufficient  range  to  permit 
changing  apparent  angle  of  arrival  by  0.  5 milliradians  (polar  and/or 
declination  components)  with  resolution  and  resetability  of  approximately 

0.  5 iLiradians.  A soldering  iron,  with  the  tip  removed  to  provide  a cavity 
six  diameters  deep,  was  used  as  a 'black'  radiance  source.  [Note  that  only 
AC  signals  are  recorded  with  low  and  high  frequency  corners  at  0.  2 Hz  and 
4 kHz,  respectively.  That  is  to  say:  a non-fibrillating  and  stationary  star 
image  produces  no  recorded  signal.  Any  voltage  changes  recorded  during 
stellar  observations  represent  object  scintillation  or  momentary  displace- 
ments from  the  mean  subaperture  image  position.  ] For  calibration  it  was 

necessary  to  modulate  the  source.  Accordingly,  a room  temperature  shutter 

-3  2 

driven  at  45  Hz  was  positioned  between  the  point  (10  cm  ) aperture  and 
the  hot  cavity.  Source  temperature  was  adjusted  to  produce  an  image  flux 
level  in  the  same  decade  as  that  anticipated  from  the  planet  Mercury.  Tape 
recordings  were  made  with  east-west  and  north-south  displacements  of  0, 

1,  2,  4,  8,  12  and  20  microradians  from  the  reference  line  of  sight  for 
AMTA  filters  #1,  3,  4,  6 and  9.  The  image  was  maintained  at  each  dis- 
placement for  one  full  minute  to  permit  precise  averaging.  Exchanging  a 
small  frosted  lamp  for  the  soldering  iron  produces  a visual  image  for  the 
AMTA  boresight  TV  camera  display  which  could  then  be  marked  with  a 
reference  fiducial  representing  the  quiescent  line  of  sight.  Planet  and 
star  images  were  brought  to  this  mark  during  the  course  of  atmospheric 
observations  to  circumvent  mount  model  and  tracking  uncertainties. 


6.  For  recent  examples  see;  G.  Ashley;  Tracking  et  al.  Trials,  Tech 
correspondence  file  nos.  A/OP-346,  A/TB-059. 
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Figures  8(a)  and  8(b)  relate  source  temperature  and  equivalent  entrance 
aperture  irradiance  to  line  voltage,  the  measured  parameter.  Figures  9(a) 
and  9(b)  relate  detector  signal  to  image  displacement  for  representative 
detectors.  Predictable,  for  small  displacements,  signal  change  is  linearly 
related  to  subaperture  image  displacement  and  detectors  are  relatively 
insensitive  to  wander  components  normal  to  the  line  joining  quadrature 
detector  and  subaperture  image  centroids.  The  voltages  measured  are 
functions  of  source  intensity  and  wavelength  as  well  as  image  displacements: 

AV  = E(X)  [S^(o)  - S^(Ar)]; 

where  E(X)  is  the  effective  in-band  entrance  aperture  irradiance  for  filter  X 
while  S^(Ar)  represents  response  in  volts  • w"  • cm  for  the  nth  detector 
with  the  subaperture  spot  displaced,  Ar,  from  its  quiescent,  (o),  position. 
Wavelength  has  two  mechanisms  for  changing  instrument  response:  (1)  the 
spectral  quantum  efficiency,  of  the  detector;  and  (2)  diffraction  which 
changes  the  spot  size.  Figure  10  illustrates  how  spot  size  is  predicted^^^ 
to  change  with  wavelength  near  focus.  S^(Ar)  includes  and  the  area 
common  to  spot  and  detector  weighted  by  intensity  contours  of  the  former 
and  sensitivity  contours  of  the  latter.  Figure  11  shows  how  normalized 
signal  should  relate  to  image  displacement  for  a uniform  detector  and  a 
shortwave  spot  (X  to  the  left  of  the  minimum  in  Figure  10).  It  is  clear  from  j 

Figure  11  and  the  calibration  data  (Figure  9(a))  that  random  spot  motions  j 

about  the  quiescent  position  will  produce  a non-zero  average  signal  history. 

The  long  period  average  value  of  the  voltage  history  is,  in  principle,  a = 

measure  of  the  average  tilt  of  the  wavefront  and,  of  course,  any  quiescent  , 

misregistration  of  the  image  illumination  centroids.  J 

The  calibration  procedure  produced  a table  of  tilt  responsivities, 

S(Ar),  for  each  detector/spectral  option  combination  against,  in  effect,  a 

standard  irradiance  source.  To  reduce  observational  data  to  tilt  structure 

size,  it  was  necessary  only  to  determine  the  in-band  aperture  irradiance  i 


7.  See  for  example:  E.  Wolf,  Proc.  Roy.  Soc.  A,  204,  p.  533  (1951). 

8.  For  example:  N.  Carleton  (editor).  Methods  of  Experimental  Physics. 
22,  Part  A,  p.  453  (1974). 
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TYPICAL  POLAR  DETECTOR  b)  TYPICAL  DECLINATION  DETECTOR 


b)  Typical  Declination  Detector 
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Figure  10  80%  Encircled  Energy 
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Figure  11  Predicted  Displacement  Response,  11  fjxn 


produced  by  the  one  star  a • Bootis,  and  the  three  planets  observed.  The 
former  comes  directly  from  astronomical  surveys.  ' ' Computation  of 
irradiance  from  solar  system  objects  is  a tedious  but  straightforward  task 
with  marked  differences  in  approach  for  each  planet.  Mercury,  for  example, 
has  no  atmosphere  and  travels  in  a rather  eccentric  orbit.  Its  subsolar 

point  has  a grey  body  temperature  of  493°K/^7R.  r.  and  an  emissivity  of 

(9)  ■*' 

0.  82;  where  y is  its  distance  from  the  Sun,  in  astronomical  units, 

at  the  time  of  the  observation.  Its  dark  side  temperature  is  taken  to  be 

111°K  from  the  measurements  of  Murdock  and  Ney.  The  illumination  • 

geometry  can  be  calculated  from  American  Ephemeris^^ ^ ^ tabulations. 

The  subsolar  point  on  Mercury  was  not  visible  during  our  observations  and 

a limb  cooling  law^^^^  T((^)  = cos^'^^  (|)  was  used  to  calculate  radiant 

intensity;  where  T(^  ) is  the  temperature  of  any  area  increment  at  a planet 

centered  angular  distance,  <^,  from  the  subsolar  point  at  temperature 

"^ssp  ~ f(date).  The  calculated  signatures  for  Mercury  during  our 

observations  is  shown  in  Figure  12. 

Venus,  with  its  dense  atmosphere,  is  essentially  a uniform  tempera- 
ture disc  and  its  long  wave  signature  is  a function  of  range  only  and  not 
solar  phase  angle.  Mars  has  a tenuous  atmosphere  and  represents  compli- 
cated intermediate  case  wherein  both  the  surface  and  atmosphere  contribute 
to  the  signature.  Fortunately,  owing  to  its  solar  distance.  Mars  can 
appear  only  slightly  gibbous  from  Earth  and  the  planet  is,  moreover,  a 
very  popular  subject  for  infrared  observations.  Thus,  its  illumination 
geometry  does  not  alter  its  observed  signature  greatly  and  many  measure- 
ments of  its  irradiance  are  available  in  the  literature. 


9.  O.  L..  Hansen,  Ap.  J.  190,  p.  715  (June  1974). 

10.  T.  L.  Murdock,  E.  P.  Ney:  Science,  170,  p.  535  (October  1970). 

11.  USN  Obs.  The  American  Ephemeris  and  Nautical  Almanac  for  1976. 
U.  S.  G.  P.  O.  (1976). 

12.  Mercury  is  similar  to  our  moon;  see  Saari  and  Shorthill  in:  AAS  Sci. 
and  Technology  Series.  13.  p.  57  (1967). 
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The  LOWTRAN'  ' Program  was  used  to  calculate  atmospheric 
transmission  and  to  trim  effective  center  wavelength  and  bandpass  for 
AMT  A spectral  options.  Passband,  AX,  and  center  wavelength,  ^g££»  are 
given  by 


R(X) 


AX  = 


max 


f 


S(X)W(X,  T)t  (X)  dX 
a 


(1) 


R = 0.  005  points 


and 


eff 


A/ 


X R(X)  dX 


(2) 


respectively,  where  S(X)  is  the  throughput  and  response  of  optical/sensor 

components  combined  and  given  by  manufacturer's  measurements;  T (X) 

a 

is  atmospheric  transmission  given  by  LOWTRAN  3;  and,  W(X,  T)  is  the 
Planck  radiance  function  for  source  temperature  T°K.  While  source 
temperature  generally  does  not  impact  AMTA  spectral  parameters  strongly, 
in  this  instant  it  did:  Filter  #9  has  a passband  1.  1 /Ltm  wide  centered  at 
3.  8 fjxn  for  the  planet  Mercury  and  shifts  to  AX  = 0.  6 fjxn  centered  ^g££  =1.6  /jm 
when  a • Bootis  (4900°K)  is  observed. 

A set  of  sicteen  detector  angular  responsivity  functions,  was 

derived  from  the  calibration  data  for  each  of  the  five  spectral  options. 


" <R>  £ = 


E (X)  F mS 

2JB 

r(X)  AF 


6 g 

6 V 


microradian  • watt  • cm" 


volt"^  (3) 


Sixteen  detectors  encompass  the  four  subaperture  images.  ”^R>r  is  angular 

th  th  ^ 

responsivity  for  the  n detector  and  f filter  option.  r(X)  if  the  image 

radius;  F is  the  telescope  focal  length  (960.  1 inches,  for  this  experiment); 


13.  See  for  example:  J.E.A.  Selby  and  R.  A.  McClatchey,  Atmospheric 
Transmittance  from  0.  25  to  28.  5 ijjn  LOWTRAN-3,  AFCRL-TR-75- 
0255  (May  1975). 


29 


m is  magnification  of  the  reimaging  optic,  while  S and  AF  are  scaling 
distances  shown  in  Figure  13.  E^p(X)  is  the  entrance  aperture  equivalent 
inband  irradiance  from  the  test  source.  To  scale  the  cross  correlation 
between  i^  and  detectors  during  planet  Mercury  observations,  for 
example,  the  convolved  voltage  histories  are  multiplied  by 

2 XI 

Ej^p  to  yield  squared  microradians  of  tilt. 

Scales  for  the  scintillation  or  'twinkle'  autocorrelation  functions 
were  established  by  squaring  the  sum  of  peak  voltages  put  out  by  four 
detectors,  comprising  a quadrature  set,  with  the  reference  source  on-axis. 
The  sum  signal  is  relatively  insensitive  to  image  position.  This  is  illus- 
trated in  Figure  14  which  presents  calibration  data.  The  lower  (staircase) 
trace  results  from  changing  common  area  between  the  illuminated  sub- 
aperture projection  and  the  detector  pair  shown,  as  the  point  image  source 
is  translated.  The  upper  trace  records  the  output  of  all  four  detectors  at 
the  same  time.  Signal  changes  at  extreme  excursions  are  due  to  convolution 
of  spot  intensity  and  differing  detector  sensitivity  contours.  The  extreme 
excursion  moves  all  spots  through  approximately  half  their  individual  radii. 
The  data  shown  is  for  spectral  option  #9  P which  is  centered  at  3.  6 /jtn. 
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Figure  14  Displacement  Calibration 
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3.0  DATA  DISCUSSION 


The  atmospheric  and  calibration  data  obtained  between  8 and  16 
Spetember  1976  was  returned  to  the  Everett  Laboratory  for  analysis.  Mon- 
itoring equipment  used  during  observations  gave  some  indication  the  signals 
recorded  in  the  1.7,  10.  5 and  11.5  micron  bands  were  correlated  between 


certain  detectors. 


The  existence  of  correlation  was  estimated  in  a 


rather  crude  way;  The  voltage  of  four  detectors,  monitoring  different  sub- 
apertures, were  written  out  simultaneously  with  a multi. pen  chart  record. 
The  traces  appear  to  track  more  frequently  when  a target  is  in  the  field  than 
when  we  'stepped -off  to  look  at  background  (see  Figure  15);  thus  it  wasn't 
a totally  blind  experiment.  At  Everett  the  signals  from  pairs  of  detectors 
were  band  limited  and  correlated  with  results  shown  in  accompanying  illus- 


trations . 


Figure  16  is  a plot  of  the  cross  correlation  function; 


R.R. 


r(T)  = 


• K..  r 


V.{t  - t)  dt,  where  -0.  5 sec  < T ^ 0.  5.  (4) 


The  integration  time,  T,  varies  between  two  and  just  over  four  minuteb 

for  this  and  the  plots  which  follow.  Scale  factors  R.  and  R.,  for  the  i*^” 

.th  ^ ^ 

and  j detector  were  determined  with  a translating  point  source  as  des- 
cribed in  Section  2.  0.  The  heavy  trace  results  from  observations  of  planet 
Venus  in  the  10.  5 micron  band  (filter  #5)  at  approximately  47°  elevation. 

The  light  trace  results  from  background  data  obtained,  a few  nunutes  lat^r, 
by  shifting  the  telescope  line-of-sight  eighty  arcseconds  north.  A cartoon 
in  the  figure  shows  the  detectors  whose  signals  were  correlated  and  the  fact 
that  they  share  a common  subaperture.  The  negative  correlation  observed 
near  t = 0 must  be  due  to  east-west  motion  of  the  light  spot;  that  is  to  say, 

the  polar  component  of  angle -of-arrival  change.  The  peak  variance  of  the 

2 

tilt  amounts  to  approximately  3.2  j^^radians  . This  measurement  was  obtained 


(14)  R&D  Status  Report,  27  October  1976,  Contract  F30602-75-C-0235. 
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with  an  eight-inch  (8ub)aperture.  We  attribute  the  background  trace  to 
microphonic  noise  (refrigerator  thumping)  and  large  scale  radiant  atmo- 
spheric structure.  The  notch  in  the  heavy  trace,  at  t = 0 is  probably  also 
microphonic  and  perhaps  intensity  scintillation  which  will  be  discussed 
presently.  Figure  17  shows  what  the  orthogonal  detectors  observed  during 
the  same  time  interval.  The  cyclic  pattern  in  the  target  trace,  for  T 
> 0.  06  secs,  is  due  to  a dysfunctioning  declination  axis  servo  drive.  The 
mount  has  a tendency  to  oscillate  with  a frequency  of  five  to  ten  hertz. 

Wind  buffeting  excites  the  critical  resonance  as  shown  in  Figure  18.  On 
14  September  the  instruments/balance  weights  configuration  and  the  servo 
adjustment  produced  a tendency  to  oscillate  at  eight  hertz. 

Figures  19  and  20  are  plots  of  correlation  between  signal  histories 
of  detectors  monitoring  different  subapertures.  The  heavy-trace  data  leads 
at  once  to  the  notion  that  relatively  long-lived  refractive  structure  is  flow- 
ing across  the  aperture  at  a fix  rate.  Negative  delay,  T,  for  r(T) 
means  tilt  events  propagate  from  east  to  west.  The  apparent  wind  velocity 
and  direction  may  be  calculated  from  this  data  because  the  detector  array 
geometry  and  its  projection  into  telescope  entrance  space  is  accurately 
known. 

For  example,  referring  to  the  figure 
at  the  left,  wherein  dj^  and  are 
N distances  between  detector  projections 

— and  vT  is  the  wind  vector,  it  is  clear 

\ by  inspection  that: 


dj  cos  (j) 
A t. 


d^  cos  (Aa  + <j>) 


from  which  is  foUo'vs  that: 


(p  = arc  tan 


- ^)] 


The  Table  2 shows  wind  speed  and  direction  calculated  from  11.5  micron 
data  obtained  9 September. 
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TABLE  2.  FILTER  #4,  11.5  /jm,  9 SEPTEMBER  1976 


I 


t 


I 

I 

* 

] 

i 


I 


Wind 

Data  Run  # Detectors  Velocity  Direction 


92 

21 

X 

6 

8.  1 m/sec 

90 

21 

X 

5 

6.  86 

89 

88 

21 

X 

7 

6.  34  I 

1 116.19° 

21 

X 

25 

6.  34  ! 

87 

21 

X 

17 

6.  69 

86 

11 

X 

17 

(5.44) 

84 

21 

X 

24 

7.  21  1 

83 

21 

X 

9 

(9.  50)  [ 

116. 05° 

82 

21 

X 

11 

7.05  Ij 

1 

116.010 

35 

21 

X 

24 

7.  21  ' 

34 

21 

X 

5 

6.  59 

31 

19 

X 

16 

7.  53 

Average 

(10  pts.  ) 

6.  99  + 0.  17  meters/second 

Calculated  angles  and  velocities  are  the  result  of  four  minute  observations. 
The  experimenter's  notes  indicate  local  winds  were  '15  mph  from  the  south- 
east' at  this  time.  As  expected,  the  air  eddies  producing  the  refractive  dis- 
turbances are  evidently  close  to  the  telescope.  NOAA  balloon  survey  records, 
taken  l-l/Z  hours  later  above  Hilo  (Hawaii)  100  miles  southeast,  show  transi- 
tion between  surface  trade  winds  (general  northeasterlies)  and  the  jet  stream 
from  the  west  is  already  well  developed  at  5 km  altitude.  AMOS  Observatory 
is  at  3.  05  km  altitude,  NOAA  data  from  Lihue  (Kauai)  250  miles  northwest, 
indicates  the  transition  occurs  at  6.7  km. 

The  11.5  micron  data  was  obtained  by  observing  the  planet  Mercury 
at  approximately  sixty  degrees  elevation  and  22:30  hours,  9 September  U.T. 

Figures  21,  22,  23,  24,  25  and  26  show  correlation  computed  for  the 
detector  pairs  illustrated  in  each  figure.  We  can  see  wind  refractive  struc- 
ture propagating  across  the  line-of-sight  and  the  following  features  are  evi- 
dent: 


* 
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a)  Structure  is  of  the  order  of  half  the  size  of  the  telescope. 

b)  The  refraction  features  may  be  highly  elliptical,  perhaps  even 
striated  like  water  waves. 

We  have  insufficient  data  to  establish  whether  or  not  this  is  typical. 

Referring  to  Figures  21  and  22  note  from  the  heavy  trace  that  maxi- 
mum correlation  is  displaced  from  T - 0 even  though  the  detectors  share 
a common  subaperture;  from  this  we  deduce  there  is  appreciable  wsivefront 
curvature  over  20  cm  of  aperture.  The  full  width  of  the  negative  maxima 
are  of  the  order  of  0.  11  seconds  which,  given  the  wind  velocity,  leads  to 
the  notion  the  refractive  ripples  have  a crest-to-crest  length  of  the  order 
of  80  cm.  Figures  27  and  28,  the  transforms  of  the  time-domain  plots, 
show  the  spectral  density  of  these  correlations.  Tilt  variance  density  is 
plotted  against  spatial  frequency  by  using  the  cosine  fall  of  the  calculated 
wind  vector  along  the  line  joining  the  appropriate  detectors.  The  ellipticity 
is  clearly  evident.  At  high  spatial  frequencies  the  spectral  density  is  pre- 
dicted to  diminish  with  eight-thirds  power  of  frequency.  The  data  fits  theory 
fairly  well  in  this  respect. 

More  striking  evidence  for  structure  size  is  shown  in  Figures  29 
and  30  which  present  10.  5 micron  data  obtained  on  14  September.  During 
this  observation  the  apparent  wind  speed  and  direction  calculated  as  just 
described  yields  3.  16  meters  per  second  from  a heading  of  128°;  (this 
corresponds  nicely  with  experimenter's  notes:  '5-7  mph  SE  wind').  Both 
plots  shown  distinct  maxima  ranging  0.  09  to  0.  35  second  apart:  the  cor- 
responding crest-to-crest  distance  appears  to  be  0.28  to  1.  1 meters.  Each 
plot  is  the  result  of  averaging  131  contiguous  one  second  cross-correlations. 

A qualitative  verification  of  the  128°  wind  heading  is  at  once  evident 
in  Figures  31  and  32  wherein  the  traces  are  the  result  of  averaging  500  half- 
second  correlations  of  Venus  data  obtained  with  the  detectors  shown.  Wind 
direction  is  given  in  each  figure  and  the  sign  of  maximum  correlation  is 
readily  predicted;  the  wind  quarter  can  also  be  deduced  from  Figure  31  be- 
cause occurs  at  T < 0.  Negative  t means  the  signal  history  of 

the  corner  detector  must  be  delayed  to  produce  a match. 

The  reproducibility  of  these  results  may  be  judged  from  Figure  33 


which  illustrates  two  data  runs  separated  by  7"°40^.  Each  trace  is  the 
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average  of  256  one-second  correlations,  between  10.5  micron-band  Venus  ! 

signal  histories  for  the  cross-hatched  detector  pair.  The  crest-to-crest  ; 

measurement  is  reasonably  reproducible  at  1. 2 meters  but  the  width  of  the  < 

maximum  changes  from  40  to  50  cm. 

The  results  with  filter  #3  (8.7  microns)  show  an  unexpected  depar- 
ture from  10.5  and  11.7  micron  data;  appreciably  less  refraction  effects 

(Hartmann  spot  wander)  are  observed  as  illustrated  in  Figures  34  and  35 

-6  2 

The  peak  tilt  variance  is  less  than  0.  5 x 10~  radians  only  ten  minutes 
before  the  11.5  micron  observation  (Figure  21).  We  have  no  explanation 
for  this. 

Figures  36  and  37  illustrate  3.8  micron  computational  results  ob-  , 

tained  from  54°  elevation  Mercury  data  on  9 September.  Correlated  tilt 
diminishes  rapidly  with  distance  in  aperture  space. 

Some  observations  were  made  with  filter  #9  to  obtain  a short  wave 
data  point.  (Filter  #9  is  simply  fused  quartz;  it  was  included  in  the  AMTA 
color  wheel  to  take  advantage  of  the  intrinsic  peak  for  germanium  detectors 
which  extends  from  the  visible  to  approximately  1.8  microns  and  to  permit 
diagnostic  testing  with  visible  stars.]  Arcturus  (a  • Bootis)  at  80°  elevation 
was  used  as  a light  source.  When  the  color  temperature  of  this  star  is 
convolved  with  the  spectral  transmission  of  atmosphere  and  instrument 
response,  a band  approximately  0.6  microns  wide  centered  at  1.6  microns 
results.  Tilt  correlation  for  1.6  microns,  calculated  from  data  obtained 
15  September,  is  shown  in  Figures  38,  39  and  40.  Calculating  apparent 
wind  speed  as  described  earlier  yields  the  correlation  scale  (in  space  in- 
stead of  time)  shown  in  Figure  40. 

Data  obtained  at  4.8  microns  and  18.8  microns  requires  special 
treatment  and  wasn't  reduced.  Figure  41  assembles  results  obtained  thus 
far  onto  a single  page.  The  plots  illustrate  how  wavefront  tilt  correlation 
between  detectors  scales  with  detector  separation  and  wavelength.  The 
straight  lines  represent  least  sqviares  power  law  fits  to  points  computed 
from  the  data.  There  appears  to  be  some  trending  with  wavelength. 

As  explained  in  Section  2.0,  amplitude  scintillation  effects  are  em- 
bedded in  the  data  and  are  expected  to  produce  positive  signed  correlation 
between  temporally  aligned  detector  signal  histories.  By  summing  together 
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Figure  37  Filter  #1  Data,  Rxjn  56 
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Figure  38  Filter  #9,  Data,  Run  97 
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Figure  41 


Tilt  Scaling  Summary 


the  four  detectors  of  a quadrature  set  we  obtain,  in  principle,  the  total 
intensity  history  for  a single  subaperture.  Owing  to  differences  in  the  sen- 
sitivity contours  of  independent  detectors  this  isn't  precisely  true  as  illus- 
trated in  Figure  14.  The  figure  shows  calibration  data  which  results  from 
translating,  in  steps,  a point  radiance  source  in  the  Cassegrain  image  sur- 
face. Translating  the  source  produces  motion  of  the  subaperture  pattern 
projected  onto  the  array  and,  consequently,  the  signal  changes  observed  in 
the  lower  (staircase)  trace.  This  trace  is  the  sum  of  signals  from  the  two 
detectors  shown  cross-hatched  in  the  figure.  The  upper  trace  is  the  sum 
of  four  detectors,  made  for  the  same  calibration  run,  and  demonstrates  the 
expected  insensitivity  to  spot  motion  except  for  extreme  excursions.  The 
abscissa  scale  is  in  j[i radians.  Since  excursions  in  excess  of  12  pradians 
were  infrequently  encountered  during  the  observations,  variations  recorded 
for  the  four  set  are  probably  not  due  to  close-in  refractive  structure.  Sam- 
ple planet  Mercury  variations  are  shown  in  Figures  42  and  43  which  are 
plots  of  the  autocorrelation  function  for  the  sum  of  the  four  detectors  illus- 
trated. Ordinate  scales  are  normalized  to  total  inband  irradiance  and  the 
traces  are  the  result  of  4 17  signal  averages.  The  peaJc  'twinkle'  ranges 

between  2 and  3%  for  these  10.  5 jjm  and  11.5  jjm  Mercury  observations. 

The  spectral  density  distribution  is  given  by  the  transforms  plotted  in  Fig- 
ures 44  and  45. 
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APPENDIX  A 


ATMOSPHERIC  AE'^ORPTION 


1 . 0 SCOPE  AND  BACKGROUND 


This  addendum  describes  work  performed  and  some  interim  results 
achieved  through  February  1976  on  the  IR  Atmospheric  Measurements  Pro- 
gram. 

The  program  consists  of  a series  of  instrumented  observations  in- 
tended: a)  to  characterize  certain  optical  properties  of  earth's  atmosphere 
and  b)  to  demonstrate  the  feasibility  of  routinely  determining  the  thermal 
signatures  of  deep-space  earth  satellites  with  an  existing AMOS. 
Work  was  begun  by  AERL  in  July  1975  in  accordance  wii^j^^^^^^^tract 
F30602-76-C-0235  as  directed  by  DARPA's  agent:  RADC^^ffll^^^^B^xli- 
fied  program,  was  scheduled  to  be  completed  May  1977. 


(The  reverse  of  this  page  is  blank) 
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2.  0 LONG  INTEGRATION  SIGNAL  PROCESSING 

Modifications  to  the  signal  conditioning  and  recording  facilities  of  an 

. . n 2) 

existing  long  wave  infrared  sensor,  AMTA,  ’ were  aiithorized  with  the 
intent  of  enabling  observations  of  faint  deep-space  satellites  to  be  made  with 
it  routinely.  AMTA,  like  all  LWIR  point -radiometers  used  under  atmosphere 
for  space  observations,  is  limited  in  precision  by  fluctuating  emission  from 
both  warm  atmosphere  and  optics.  At  AMOS  these  two  sources  produce  four 
to  five  decades  greater  flux  level,  at  the  detector  surface,  than  typical  tar- 
gets. The  line  of  sight  of  one  or  more  AMTA  detectors  is  deviated  by  an 
oscillating  mirror  to  alternately  encompass  and  exclude  the  target.  AMTA's 
scan  mirror  is  driven  to  follow  a 45  Hz  square -wave  which  toggles  the  line 
of  sight  24  arc  seconds,  an  angular  distance  corresponding  to  the  detector 
spacing.  In  the  absence  of  a target,  thermal  eddies  in  the  atmosphere  affect 
a noise -like  detector  output  which  must  be  integrated  to  produce  a smooth 
mean  value  so  the  target,  when  present,  will  obtrude. 

Although  the  spatial  distribution  of  radiant  structure  in  the  atmosphere 

was  of  intense  interest  to  the  IR  community  in  the  late  1950's,  rather  little  is 

(3) 

available  by  way  of  study  reports  today.  The  Weiner  spectrum'  is  one  en- 
during description;  it  predicts  the  power  spectral  density  of  noise,  produced 

by  a slot  scanning  the  atmosphere,  will  vary  inversely  with  the  square  of 

(4) 

spatial  frequency.  It  has  been  suggested'  ' that  a corner  ought  to  occur  be- 
tween one  and  ten  cycles  per  radian,  below  which  the  spectrum  is  uniform. 
This  description  suggests  that  the  noise  recorded  during  an  observation  ought 
to  be  highly  correlated  between  detectors  which  are  separated  by,  at  most, 

550  ^.radians.  The  background  is  expected  to  trend  monotonically  across  the 
array  and  first  order  differencing  between  adjacent  'identical'  detector  chan- 
nels should  eliminate  the  influence  of  atmosphere. 


3. 

4. 


Ibis  reference  will  be  made  available  to  qualified  military  and  government 
agencies  on  request  from  RADC  (OCSE)  Griff iss  AFB  NY  13441. 

This  reference  will  be  made  available  to  qualified  military  and  government 
agencies  on  request  from  RADC  (OCSE)  Griff iss  AFB  NY  13441. 

University  of  Michigan  report  #IRIA  2389 -7 -X  (July  1957). 

D,  Korff,  private  communication. 


Figure  A-1  shows  how  this  is  accomplished  with  AMTA,  Owing  to 
the  toggling  motion  of  the  scan  mirror  and  the  large  scale  structure  of  back- 
ground granularity,  most  of  the  detectors  (for  this  illustration)  output  45  Hz 
signals  whose  components  are  in  phase  between  detectors.  Detector  channel 
#7  is  subtracted  from  channel  #13,  therefore,  background  contribution  from 
the  two  interfere  destructively;  if  the  radiant  gradient  is  uniform  and  the 
detectors  are  identical,  cancellation  is  complete.  On  the  other  hand,  target 
signals  from  the  two  are  out  of  phase  to  begin  with  so  subtracting  #13  - #7 
sums  target  contributions.  Fine  atmospheric  strvicture,  photon  arrival  fluc- 
tuations, amplifier  and  detector  self-noise  is  uncorrelated  between  channels; 
therefore  using  a pair  of  detectors  improves  S/N  by  . Unless  the  detec- 
tor channels  have  identical  responsivities,  there  will  be  a non-zero  average 
for  the  difference,  #13  - #7,  even  in  the  absence  of  a target.  To  overcome 
this  effect  the  telescope  line  of  sight  is  shifted  to  put  the  target  on  a second 
pair  of  toggling  detectors  for  an  equal  amount  of  observation  time.  The  change 
in  the  mean  value  of  the  signal  history  for  each  pair,  produced  by  moving  the 
LOS,  is  attributed  to  target  and  no  observation  time  is  lost  while  background 
is  measured  because  the  target  is  always  in  view. 

Satellites  and  R/V's  at  slant  ranges  of  2500  km  or  less  can  be  ther- 
mally characterized  with  a few  seconds  of  data  while  small  objects  in  syn- 

3 

chronous  orbit  will  require  10  seconds  or  longer.  Satellites  in  Keplerian 
orbits  provide  the  required  observation  time.  (In  principle,  synchronous 
satellites  can  be  observed  indefinitely.  ) Unfortunately,  AMOS  experience 
teaches  that  background  statistics  are  apt  to  be  non- stationary  over  periods 
longer  than  an  ho^ir.  Double  differencing  using  'tare'  detectors  helps  to 
compensate  for  long  term  trending  of  the  mean  gradients  of  atmospheric 
radiance. 

Figure  A-2  shows  the  organization  of  signal  conditioning  and  recording 
hardware.  Four  such  two-detector  signal  processors  comprise  the  long  term 
averaging  system.  The  band  pass  filter,  illustrated  in  the  figure,  is  optional. 
The  analog  to  digital  converter  is  a precision  unit  with  l6-bit  resolution  which 
reads  the  first  integrator  45  times  per  second.  Two  types  of  software  are 
available:  FAST  and  SLOW  modes.  For  the  first  named,  the  summation 
register  acts  as  a buffer  storage  only,  and  45  readings  for  each  of  four  pro- 
cessors is  recorded  on  nine-track  digital  tape  each  second.  The  SLOW  mode 


permits  real-time  TTY  interaction  and  recording  with  some  sacrifice  in  time 
resolution:  The  average  of  45  readings  is  recorded  on  tape  once  each  second 
for  each  processor.  The  TTY  will  printout  the  average  of  the  difference  be- 
tween (selected)  pairs  of  processors  once  each  (selected)  averaging  time. 

The  minimum  averaging  time  for  this  mode  is  ten  seconds;  the  maximum 

4 

is  10  seconds. 

During  the  course  of  long  observations  the  background  signal,  out- 
putted by  the  first  integrator,  has  changed  by  as  much  as  a volt  per  hour 
-20  -2 

(trending)  while  10  w • cm  target  irradiance  (a  measuring  goal)  pro- 

_5 

duces  1.7  to  2.0  x 10  volts.  Seventeen  bits  are  required  to  encompass 
this  range.  The  l6-bit  converters  thus  provide  a 2 x 10  w • cm  un- 
certainly per  reading.  In  principle,  the  mean  value  of  a set  of  readings  may 
be  determined  with  greater  than  the  encoding  resolution  by  averaging.  For 
example:  if  two  value -adjacent  (separated  by  one  least  significant  bit)  digital 
estimates  exactly  straddle  the  actual  mean  value  the  least  significant  bit  will 
be  'lit'  for  exactly  50%  of  the  sample  readings.  The  average  of  sixteen 
readings  should  yield  one  additional  bit  resolution  with  90%  confidence. 

Sixteen-bit  precision  converters  are  available  from  two  manufacturers.  The 
only  available  converter  with  higher  resolution  is  a state-of-the-art  oven 
stabilized  digital  voltmeter,  manufactured  by  Hewlett  Packard,  which  offers 
one  microvolt  resolution  over  a 1 volt  range  (equivalent  to  21  bits).  One 
such  instrument  was  incorporated  into  the  system  to  monitor  the  performance 
of  the  other  converters.  Its  BCD  output,  which  represents  a 1.  1 second  av- 
eraged reading  is  recorded  once  every  two  seconds  along  with  the  other  data. 

The  HP  voltmeter  is  generally  left  in  parallel  with  the  1 6-bit  converter  moni- 
toring the  detector  pair  closest  to  the  nominal  telescope  line  of  sight. 

Data  processing  includes  calculating  standard  deviation,  mean  value 
and  least-squares  fit  to  first  order  trending  with  time  for  the  average  values 
obtained  per  telescope  LOS  dwell  time.  This  is  typically  ten  seconds  so  the 
16-bit  value  is  the  average  of  450  readings  while  the  HP  value  is  the  (second) 
average  of  only  five  readings.  Nevertheless,  the  HP  values  generally  have 
10%  smaller  normalized  mean  deviation  for  signal  flux  near  threshold  for 
five  minute  observations. 
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No  manufacturer  will  describe  the  performance  of  an  encoder,  a 
voltmeter,  or  A/D  converter  beyond  its  LSB  resolution.  We  tested  the  16- 
bit  units  with  signals  incremented  in  steps  smaller  than  LSB  during  system 
tests  and  again  during  operational  tests  at  AMOS.  The  units  did  exhibit  the 
ability  to  interpolate. 

' The  key  limitation  to  long  term  averaging  is  the  accurate  identifica- 

tion of  all  sources  of  trending.  For  this  reason  the  detector  bias  voltage  bus 
and  the  'absolute'  value  of  detector  signals  is  impressed  on  the  same  nine- 
track  tape.  Changes  in  apparent  sky  brightness  can  be  deduced  from  the 
absolute  value  history.  Effective  quantum  gain  for  the  photoconductors  is  a 
function  of  detector  bias  which  controls  charge  mobility.  Although  the  bias 
bus  is  reasonably  well  regulated,  changes  have  been  observed  during  long 
integration  missions. 

The  nine -track  record  also  includes  radiometer  filter  number,  Uni- 
versal time  from  the  observatory  time  code  generator,  three  event  flags  and 
a typed-in  initializing  statement  which  identifies  the  mission  name,  date  and 
detector  channel  assignments. 

Calibration  is  accomplished  by  directing  exitance  from  a chopper 
modulated  reference  black  body  into  the  telescope  entrance  aperture.  The 
chopper  is  phase  locked  to  the  system  46  Hz  reference  signal  which,  at 
other  times,  drives  the  radiometer  oscillating  mirror.  The  reference  source 
is  mounted  on  the  secondary  mirror  support  structure  of  the  1. 2 meter  tele- 
scope. A remotely  commanded  mirror  directs  black  body  energy  into  the 
telescope  as  required.  In  practice,  the  dwell  time  for  this  mirror  is  adjusted 
to  correspond  to  the  dwell  time  employed  for  the  on  and  off-axis  detector  chan- 
nels by  the  telescope  LOS  during  the  mission.  Atmospheric  extinction  is  det- 
ermined by  observing  standard  stars  and  by  best-  ' •'ing  an  absorptance  model 
to  sky  radiance  gradient  data  obtained  in  near  real  time. 

2.  1 RESULTS 

LWIR  flux  from  certain  stars,  two  asteroids  and  a number  of  deep 
space  satellites  were  measured  with  the  new  system  in  January  1976.  Some 
results  are  presented  here: 
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Long  Integration 


Intensity 


ADC  # 

Date 

Time 

Name 

Range 

Elev. 

Filter 

w • 

ster'l 

6356 

24  January 

10:40 

Moly  1-23 

37500 

39.  5° 

5 

530 

+ 120 

7625 

23  January 

08:15 

COS  #706 

40440 

52.  5° 

5 

360 

+ 160 

6356 

22  January 

09:30 

Moly  1-23 

29900 

39.  8° 

5 

730 

+ 250 

83564 

22  January 

09:59 

DSP 

39223 

23.  0° 

5 

450 

+ 360 

83564 

28  January 

08:10 

DSP 

39270 

23.  1° 

5 

640 

+ 136 

83564 

28  January 

08:30 

DSP 

39269 

23.  1° 

6 

< 410 

The  intensities  recorded  were  obtained  with  S 600  seconds  integration. 
DSP  is  a U.S.  satellite  in  synchronous  orbit  while  the  others  are  USSR  objects 
in  highly  elliptical  half-day  orbits.  The  20  fjxn  (filter  #6)  bounding  value  listed 
for  DSP  represents  the  standard  deviation  about  the  background.  The  signal 
history  showed  no  repeatable  increase  when  the  FOV  encompassed  the  target 
during  the  observation.  The  results  of  an  earlier  observation  of  Molniya  1-19 
is  shown  in  Figure  A- 3. 

Infrared  observations  were  made  of  the  following  astronoinical  objects: 


Star  Name 
(3  And 

a Tan 

R Leo 

R Hya 

MNL  Cyg 

L Per 

6 Tan 

Y Hya 

p Her 

T Pup 


AMTA  Filter  # Objectives 

2,  5,  6 Calibration  and  extinction  determination 

2,  5,  6 " " ” 

2.  5,  6 '■  " ” 

2,  5,  6 ” " " 

5,  6 " " 

5,  6 Long  average  observation 

5,  6 " " " 

g II  II  II 

g II  tl  It 

5 " " " 
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8-13.3  RADIANT  INTENSITY 


13=0000  13=05=  13=00=  13  = 15=  13=20=  I 

UNIVERSAL  TIME 

G6986 

Figure  A- 3 Molniya  1-19  Thermal  Signature 
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Asteroids 


AMTA  Filter  # 


Objectives 


I 


Ceres 

2,  5.6 

Temperature  measurement  test 

Juno 

5.6 

M tt  M 

Vesta 

2,  5.6 

It  It  It 

The  first  five  are  bright,  familiar,  infrared  stars  which,  with  the 
exception  of  p And,  are  used  for  AMTA  calibration  in  support  of  observations 
generally  not  requiring  integration  beyond  one  second.  The  signatures  ob- 
tained for  Q Tau,  MNL  Cyg,  R Leo  and  R Hya  were  repeatable  night  after 
night  within  + 2%  and  correspond  to  catalog  values^^’  within  + 9%  for 
filters  #5  and  #2.  Correspondency  at  20  (im  (filter  #6)  was  considerably 
poorer:  + 58%  or  0.  5 stellar  magnitudes.  This  is  probably  due  to  incorrect 
atmospheric  extinction  assumptions.  Considerably  fewer  itiq  observations 
are  published  than  either  m^  or  m^.  AMOS  brightness  results  for  p And 

IN  M 

compare  well  with  recent  observations  by  F.  J.  Low  and  G.  Rieke: 


”^M 

mQ 

Te„(N/Q) 

— MfitW 

AMOS 

- 1.87 

- 2,  13 

- 2.41 

1300°K 

2740OK 

Ref  7 

- 1.97 

- 2.  05 

- 2.20 

2130°K 

5500°K 

eff 

5 jLtm 

10.  6 

21  /im 

The  temperatures  listed  are  calculated  by  fitting  the  indicated  two  color  spec- 

tral  irradiances  to  a black  body,  p And  is  listed'  ^ as  spectral  type  MO  III 

star,  which  suggests^*^^  a visual  color  temperature  of  3200°K,  and  as  a sus- 
/Q) 

pected  variable,  ' ' which  might  easily  account  for  the  different  apparent  tem- 
peratures observed. 

5.  J,  Heath,  AMOS  Star  Library.  Lockheed  Missiles  and  Space  Company, 
Inc,  (April  1974). 

6.  F.  J,  Low,  Sky  Survey.  Semi-Annvial  Technical  Report.  AFCRL  #70-0179, 
Univ.  of  Arizona  (15  March  1970). 

7.  N.  Carleton  (editor).  Methods  of  Experimental  Physics.  Vol.  12, 
Astrophysics,  p,  453,  Academic  Press,  N.Y.  (1974), 

8.  D.  Hoffleit,  Catalog  of  Bright  Stars.  3rd  Ed.  Yale  University  Observ. 
(1964). 

<♦  C.  W.  Allen,  Astrophysical  Quantities.  3rd  Ed.  p,  206,  Univ.  of  London, 
Athlone  Press  (1973). 
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The  next  five  stare  listed  in  the  observation  table  are  much  fainter. 
These  were  selected  to  demonstrate  that  familiar  visual  stars  which  are 
known  to  be  stable  and  free  of  spectral  trauma  (principally  class  V,  sun-like 
stars)  should  have  predictable  signatures  in  the  LWIR.  If  this  proves  to  be 
true,  such  stare  provide  useful  calibration  standards  now  that  long  integration 
has  produced  the  required  sensitivity.  Excepting  6 Tau,  none  of  those  ob- 
served are  listed  in  our  rather  extensive  LWIR  star  catalogs.  Reference  6, 
which  describes  observations  made  before  March  1970,  reports  m^^  = + 2.  55 


for  this  star.  Our  own  measurements  indicate  m...  = + 1.  50  and  these  results 


N 


repeat  within  0.  3 stellar  magnitudes. 

Of  the  three  asteroids  observed  two,  Ceres  and  Juno,  yield  temper- 
ature and  diameter  estimates  while  Vesta  exceeded  background  only  for 
5 jLtm  (filter  #2).  I suspect  we  were  actually  looking  at  an  my  = +8.9,  KO 
star  (GC  818).  The  apparent  surface  temperature  of  the  asteroids  Ceres 
and  Juno  was  measured  using  the  newly  installed  system.  Multicolor  data 
from  20  and  22  January  observations  were  also  used  to  determine  size  and 
surface  albedo  for  these  objects.  The  newly  developed  atmospheric  extinction 
models  were  used  for  filter  #6  data. 


radius. 

km 

AMOS 

Allen^^^ 

Range, 

Temp°K 

Ceres 

+ 7.  7 

409  + 34 

380 

3.049  x 10® 

235  + 7.4°K 

0.  86 

Juno 

9.6 

155  + 33 

100 

2.693  x 10® 

215  + 18.  30K 

0.  51 

It  is  assumed  that  these  objects  radiate  as  grey  bodies  between  8 and  22  fiv' 


with  emissivity  something  like  earth's  moon,  i.e.  : e 


~ C /XT 

In  the  regime  where  e 2'  » 1 spectral  intensity  of  thermal  radia- 

tors is  exponentially  proportional  to  inverse  temperature,  and  obversely, 
temperature  is  given  by 


0.89. 


T = 


r ^xi  1 ” ^ 

= I Kj  In  — + K2  , I^j  and  1^2  intensity  at  XI  and  X2.  (A-1) 

L X2  . 


For  AMTA  filters  #5  and  #6  Kj  = 1.8833  x lO'^  and  K2  = 5.0602  x 10~^, 


as  determined  from  baseline  measures  of  instrument  parameters,  where 
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T is  in  °K.  Spectral  intensity  (or  exoatmospheric  irradiance)  was  estab- 
lished by  instrument  response  to  known  black  body  existance  and  correcting 
for  atmospheric  extinction.  Temperatures  listed  in  the  table  are  the  aver- 
age of  four  permutations  of  two  data  sets,  for  each  color,  for  each  of  two 
detectors.  Area  is  calculated  directly  from  exoatmospheric  irradiance 
because  range  is  known.  The  objects  are  assumed  to  be  spherical. 

Albedo  is  calculated  from  the  temperature  as  follows:  for  passive 
grey  bodies  in  equilibrixim  with  solar  radiation,  temperature  goes  with  in- 
verse square  root  of  solar  distance  and  the  maximum  temperature  a grey 
body  can  reach  is  393.  6°K/n^;  where  R is  solar  distance  in  A.  U.  The 
measured  temperature,  a function  of  the  optical  properties  of  the  surface, 
is  given  by: 

A 


Qg  = a 


^TH 


(A -2) 


where  S^,  Ug,  a and  are  the  solar  constant,  effective  absorptance  of 

the  surface  for  solar  radiation,  Stefan -Boltzmann  constant  and  thermal 


(total  hemispheric)  emissivity,  respectively.  For  grey  bodies  a-  = err-rr 

o TH 


and 


(393. 6) 


max 


(A-3) 


Albedo  is  then  given  by 


1 - Ug  = 


1 - 


(393.  6) 


-TH 


{A-4) 


Observation  of  Ceres,  Juno  and  stars  for  calibrating  the  sensor 
yielded  data  for  calculating  NEFD;  some  results  are  given  here: 


NEFD,  Filter  #5 
Integration  Time 


Det.  # 

1 sec 

120  sec 

360  sec 

600  sec 

13 

18Z 

4.  9Z 

1.  IZ 

0.  77Z 

25 

13Z 

3.  2Z 

2.  5Z 

0.  76Z 

10.  American  Ephemeris  and  Nautical  Alamanac  for  1976,  U.  S.  GPO. 
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NEFD,  Filter  #6 

Integration  Time 

Det.  # 

1 sec 

120  sec 

360  sec 

600  sec 

25 

26Z 

6.  2Z 

4.  3Z 

13 

21Z 

3.  7Z 

2.4Z 

I The  noise  equivalent  fltix  density  is  determined  from  the  standard  deviation 

[ about  the  mean  for  stellar  signals.  Popular  practice  for  reporting  sensitivity 

[ entails  dividing  measured  s.  d.  by  the  square-root  of  the  number  of  readings; 

i while  this  yields  a more  satisfying  resxilt,  it  is  mathematically  correct  only 

i for  stationary  random  processes  where  the  individual  readings  are  truly  in- 

dependent.  W ith  long  term  trends  produced  by  small  local  weather  condition 
changes  one  caxmot  be  certain  the  statistical  properties  of  the  background 
^ remain  stationary.  Note  that  these  results  apply  to  single  detectors:  #13 

I and  #25,  not  combinations  of  two  adjacent  detectors  discussed  earlier.  The 

double  differencing  capability  had  not  been  tested  during  the  performance 
j period  covered  by  this  report. 
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3.  0 ATMOSPHERIC  EXTINCTION  SCALING  - TASK  4.  1 
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Thermal  observations  of  foreign  satellites  can  yield  important  mis- 
sion clues.  The  effective  surface  temperature  of  operational  satellites  may 
be  determined  by  making  multi-color  measurements  of  thermal  exitance. 

When  such  observations  are  made  from  beneath  atmosphere  extinction  cor- 
rection uncertainties  at  long  wavelengths  can  seriously  limit  the  validity  of 
conclusions.  The  principal  source  of  uncertainty  is  the  amount  of  H2O, 
which  strongly  affects  absorption  in  the  Q-window  (16-23  /jm)  region,  along 
the  line  of  sight.  On  mountain  tops,  where  abrupt  weather  modifications  are 
common,  water  vapor  content  changes  so  rapidly  that  extinction,  determined 
by  observing  stars,  may  not  be  applicable  to  a satellite  measurement  made 
an  hour  later.  During  the  course  of  LWIR  observations  we  have  recorded^^^ 
a monotonic  sky  signal  increase,  correlated  with  increasing  zenith  distanpe, 
when  the  instrument  is  scanning  clear  sky.  At  20  pm  this  signal  is  princi- 
pally due  to  thermal  emission  from  atmospheric  water  vapor  as  shown  in 
Figure  A-4.  By  measuring  radiance  gradient  an  absorption  model  may  be 
calculated.  Moreover,  if  the  radiometer  is  arranged  to  provide  am  abso- 
lute measure  of  background  radiance,  extinction  may  be  deduced  for  a satellite 
mission  even  when  there  is  no  orderly  relation  between  absorptance  and  zenith 
distance . 

To  prove  this  concept  facilities  for  measuring  and  recording  mean 
detector  resistance  for  ten  detector  channels  were  incorporated  into  the 
AMTA  system.  This  facility  shares  (see  Figure  A-5)  a tape  recorder,  rack 
space  and  a minicomputer  with  the  long  period  averaging  system  described 
earlier  in  this  report.  Resistance  (actually  d.c.  detector  output)  is  monitored, 
digitized,  multiplexed,  averaged  and  recorded  on  nine-track  digital  tape  along 
with  the  contrast  data  from  astronomical  and  satellite  observations  without 
interfering  with  or  degrading  AMTA  performance  in  any  way.  The  hardware 
was  installed  at  AMOS  in  December  1975  and  became  operational  in  January. 
Figure  A-6  is  a system  block  diagram  of  the  extinction  measuring  add-on  for 
AMTA. 
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The  d.c.  output  of  each  of  ten  detectors  is  digitized  with  12-bits  per 
ten  volts  resolution,  eight  times  each  minor  scan  cycle.  The  four  readings 
for  each  of  the  two  scan  mirror  positions  are  averaged  and  recorded  sepa- 
rately. This  produces  ninety  recorded  readings,  each  second,  for  each 
detector.  These  numbers  are  retrieved  from  tape  and  averaged  to  any  de>- 
sired  degree  of  smoothness  along  with  the  mission  data.  An  eleventh  12-bit 
A./D  converter  monitors  detector  bias  which  is  also  recorded  on  tape  ninety 
times  per  second. 

Detector  voltage  is  related  to  absolute  radiance  by  means  of  the  same 
black  body  reference  source  used  to  interpret  mission  data.  To  unwind  in- 
strument bias  and  off-set  due  to  emission  of  the  optics  the  sensor  filter  wheel 
is  rotated  to  block  all  detectors  with  an  opaque  port  refrigerated,  like  the 
filters  themselves,  to  si40°K.  The  detector  array  lies  at  the  bottom  of  a 
cylindrical  low  reflectance  light  shield  and  is  maintained  at  11°  to  14°K. 

The  array  semiconductors  will  not  respond  to  photons  less  energetic  than 
22  fjm;  the  blocked  light  shield  thus  establishes  the  'zero  photon'  detector 
resistances.  Next,  the  blocking  port  is  replaced  with  the  interference  filter (s) 
to  be  used  during  the  mission;  (the  AMTA  filter  wheel  is  remotely  commanded 
and  its  disposition  is  automatically  recorded  on  the  data  tape  once  each  second) 
and  the  telescope  dust-covers  are  closed.  The  resulting  detector  outputs 
represents  response  to  the  optical  train  exitance  with  ass\imed  unity  emissivity. 
The  calculated  emissivity  for  all  cascaded  optical  surfaces  which  are  not  re- 
frigerated is  € = 0.  205,  This  value  is  applied  to  the  difference  between  the 
dust  cover  reading  and  the  blocked  detector  reading  to  establish  what  fraction 
of  the  output,  recorded  during  the  measurement,  to  attribute  to  warm  instru- 
ment optics . 

In  January  1976  tests  were  conducted  with  the  newly  installed  hardware 
to  confirm  the  validity  of  relating  atmospheric  brightness  to  transmission. 
Measurement  procedures  were  described  in  MIOP  015  {Mission  Instruction 
and  Operational  Plan)  distributed  at  AMOS.  The  data  taken  in  January  did 
produce  a plausible  zenith  distance  related  expression  for  atmospheric  ex- 
tinction but  sufficient  'standard  star'  observations  to  fully  validate  the  model 
were  not  undertaken,  due  to  press  of  other  work.  The  model  at  20  fxm  (filter 
#6)  assumes,  since  most  of  the  water  is  in  the  first  2 km  above  AMOS,  that 
a temperature  scaled  1.4  km  above  altitude  can  be  assigned  to  the  dominant 
radiating  constituent  and  that  the  output  produced  by  looking  just  above  0°  ele- 
vation is  representative  of  total  absorption  (€  = 1.0). 
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The  data  tapes  were  reduced  in  Everett  where  software  for  computing 
extinction  from  spectral  radiance  data  was  completed  and  exercised  with  "sky 
scan"  radiance  data  obtained  from  AMOS  as  a result  of  January  activities  there. 

The  apparent  spectral  radiance  of  the  atmosphere  is  attributed  to  its 
effective  temperature  and  emissivity.  A radiance  gradient  measurement  will 
yield  emissivity  and  transmission  since  1 - e(X)  = T(X).  Transmission  as 
function  of  elevation  angle  is  given  by 


T(X)  = exp 


for  uniform  exponential  atmosphere  (A- 5) 


sin  El 


above  “ lO”  elevation  (4  air  masses  at  AMOS).  Pj^,  are  atmospheric 
pressure  at  the  altitude  of  the  observer  and  sea  level  respectively;  and  P(X) 
is  a number,  ^ 1.0,  depending  upon  the  concentration  profile  of  the  absorbing 
medium  and  whether  the  line  absorbers  are  saturated.  In  the  visible,  where 
extinction  is  attributed  primarily  to  scattering  and  the  medium  is  tenuous  for 
broad  spectral  intervals,  P(vis)  = 1.0.  In  the  LWIR  where  molecular  absorp- 
tion dominates,  and  the  line  centers  are  apt  to  be  saturated,  theory  predicts 
p(X)  ^ 0.  5.  The  radiance  data  obtained  in  January  seems  to  bear  this  out  at 
19  jim  but  not  at  10  jLltn. 

When  fully  developed  the  sky  foreground  radiance  measuring  scheme 
promises  to  permit  unwinding  the  transmission  affected  modulation,  produced, 
e.  g. , by  the  line -of-sight  crossing  patches  of  localized  water  vapor,  from 
proper  target  attitude  change  affected  signals.  Figure  A-7  illustrates  just 
such  a disorderly  sky  during  the  'Filter  #2'  data  run.  In  the  figure  detector 
output  voltage  is  related  to  zenith  distance  expressed  in  sea  level  equivalent 
air  masses.  Experiments  conducted  in  January,  however,  were  designed  to 
prove  the  general  concept  and  data  was  fitted  to  a monotonic  sky  model. 

The  data  recorded  on  magnetic  tape  22  January  1976  consists  of  a time- 
tagged  history  of  d.  c.  voltage  outputted  by  five  detector  channels  while  the 
telescope  stepped  in  elevation  from  + 4*^  to  zenith  in  accordance  with  a pre- 
arranged schedule.  Since  detector  channel  output  change  is  linearly  propor- 
tional to  radiance  change: 


OmCTOR  **15,  FILTER  •’6, 
BIAS*  -4.2S4V0LTS 
22  JAN  If  7«  ib:2o;oo  U T. 


REPORT  CLEAR  SKY 
DOME  TEMP  40*  F 
* R.K48% 


DETECTOR ’'IS,  FILTER’^? 
BIAS  • -4.2S4  VOLTS 
22  JAN  IBTS  19:20  ;OOU  T. 


REPORT  CLEAR  SKY 
DOME  TEMP  40* F 
• R.H.  4S% 


SEA  LEVEL  EQUIVALENT  AIRMASS 

-64.T4 

SEA  LEVEL  EQUIVALENT  AIRMASS 

OCTCCTOR  •!»,  riLTEA  *“9 

~ -94.76 

1 ^ iAfc 

Vw  9IA,  • -4M9  VOLTS 

1 

Q 
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JAN.  ITT*  la  St  OOUT. 

Z -W.TB 

iO 

i 94  S0 

“a  BIAS  • -4.269  VOLTS 

L 22«IAN.  1976  I9;49:00U.T. 

O 

REPORT  OCCATIONAL  \ 

< 94  M 

~ VERY  LIGHT  CLOUDS  \ 

Li 

DOME  TEMP  40*F  \/  

REPORT  CLEAR  SKY 

-94.64 

• R.H.99% 

DOME  TEMP  • 40*P 
• RH  'ASU 

-54.661 

■ 
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2.00  4.00  6.00 

SEA  LEVEL  EQUIVALENT  AIRMASS 


Figure  A-7  Air  Mass  Related  to  Detector  Output 


E(E1)  = 1 - T(X)  . 


E(E1)  = E^  1 - exp 


and  ultimately 


(A.6) 


•5^  sin  El 


In  In  - In 


Q(X)  - p(X)  In^^  sinEl^; 


which  can  be  fitted  by  least  squares  to  the  data  to  yield  a(X)  and  ^(X).  E^ 

is  detector  response  to  a black  radiator  at  the  effective  temperature  of  the 

atmosphere;  for  20  fjxn  this  is  essentially  the  voltage  for  zero  elevation.  Sky 

scans  were  made  for  filters  2,  5,  6 and  7. 

Filter  #5  (X  = 10.9  um,  AX  =4.1  um)  and  #6  (X  =18.9  uni, 
m m 

AX  = 3.3  jLim)  data  was  reduced  with  the  new  program  to  yield  the  following 

re  suits ; 

a(X)  + a p(X)  + a 

Filter  #5  0.200  + .006  0.372  + .008 

0. 610  + . 028 


Filter  #5 
Filter  #6 


Filter  #5,  r(10.9(im) 


Filter  #6,  t(18.9  /im) 


Air  Masses 


1. 109  + . on 


Zenith 

4 5°  Elev. 

0.84  + .004 

0.82  + .004 

0.41  + .004 

0.  33  + . 004 

0.  704 

0.996 

20^  Elev. 


0. 77  + . 004 


0. 18  + . 008 


2.  058 


The  values  and  scatter  for  a(X)  and  P(X)  are  the  averages  for  five 
detectors;  Eq.  (A- 5)  was  used  to  compute  transmission. 

The  sensor  was  calibrated  with  a precision  black  body  and  then  ob- 
servations were  made  of  the  apparent  10.9  fjm  and  18.9  pm  brightness  of  the 
IR  star,  NML  Cygnus,  at  47°  elevation.  Using  the  average  of  eight  exoatmos- 
pheric  brightness  values,  from  astronomical  literature,  yielded  transmission 
of  0.  77  and  0.  30  for  filters  5 and  6,  respectively;  this  should  be  compared  to 
the  45°  values  deduced  from  radiance  gradient  data.  Published  brightness 
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values  span  28%  about  the  mean.  Since  this  is  typical  for  IR  stars,  extipq- 
tion  is  more  properly  determined  by  following  a single  star  as  it  sets  through 
two  or  more  air  masses,  or  by  observing  several  "standard"  stars  at  different 
zenith  distances.  But  there  are  very  few  well  observed  IR  stars  with  N and  Q 
brightness  uncertainties  small  (compared  to  transmission  uncertainties)  listed 
in  the  catalogs.  Following  rising/setting  stars  entails  telescope  time,  which 
was  not  available  in  the  press  of  other  work,  and  in  any  event  yields  chancy  re- 
sults with  changing  weather. 

j 

Filter  #2  (5  jjjn)  and  filter  #7  (18  + 0.  8 (jim)  data  from  22  January  and 
more  data  filters  #5  and  #6  obtained  during  changing  weather  on  11  February 

, j 

1976  and  being  analyzed.  A total  of  seven  sets  of  additional  sky  scans  plus 
accompanying  calibrations  have  been  requested.  ! 

i ‘ 

■ 


Note  that  the  ratio  of  T(  1 8.  9)/t(1  0.  9)  reduced  from  radiance  measures 
and  from  star  data  agree  very  closely:  0.  39  and  0.40  respectively;  the 
apparent  temperature  of  NML  Cyg  is  better  known  than  its  brightness. 


(The  reverse  of  this  page  is  blank) 
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4.  0 MEASURING  PRECISION  CONSTRAINTS 

To  improve  immunity  to  the  AMOS  EMI  environment,  some  long 
term  drift  problems  attributed  to  cabinet  temperature  changes  and  common 
mode  response  problems,  hardware  modifications  described  during  the  19 
February  briefing  are  being  implemented.  Fabrication  is  20%  complete  and 
an  additional  twenty  man  days  of  labor  should  suffice  to  complete  this  work. 

The  sky  scan  data  produced  verification  of  the  long  reported  belief 
that  radiant  structure  in  the  atmosphere  limits  the  sensitivity  of  a fast- slewing 
radiometer.  Wind  driven  eddies  are  observed  at  AMOS  even  when  tracking 
stars;  F.  J.  Low^^^  refers  to  'sky  noise'  which  is  probably  the  same  thing. 

The  sky  scan  measuring  procedure  in  January  involved  moving  the  telescope 
up  from  the  horizon  to  zenith  in  an  orderly  series  of  stepstare  operations. 
Detector  resistance  (d.  c.  voltage)  was  recorded  for  five  channels  in  order 
to  chart  the  relation  between  apparent  brightness  and  zenith  distance  as  ex- 
plained earlier.  The  highly  amplified  a.c.  outputs  were  recorded  too.  This 
was  done  to  provide  data  relating  NEFD  to  zenith  distance,  a sort  of  no -effort 
aside  to  the  primary  measurement.  It  turned  out  that  deviation  about  the 
mean  background  was  consistently  higher  for  'step'  than  for  'stare'.  Standard 
deviation  of  the  recorded  background  was  converted  to  equivalent  entrance 
aperture  flux  by  relating  to  the  output  voltage  change  produced  by  the  black 
body  reference  source.  The  results  are  plotted  in  Figure  A-8.  The  dotted 
trace  is  a least- squares  fit  to  the  'stepping'  data;  the  solid  curve,  which  ties 
the  'stare'  points  together,  is  calculated  for  'still-air'  sky  and  shifted  upward 
a factor  of  six.  Note  that  above  40°  elevation  moving  the  telescope  produces 
a factor  of  two  more  'uncertainty',  but  conditions  wherein  motion  makes  no 
difference  do  exist  here  while,  at  lower  elevation,  motion  consistently  de- 
graded measuring  precision.  On  rare  occasions  measurements  near  zenith 
with  a fixed  telescope  have  yielded  2.  5 x Z NEFD.  The  results  shown  with 
slewing  help  to  explain  why  data  from  satellites  in  low  orbit  and  data  obtained 
at  low  elevation  is  so  unexpectedly  'ragged'.  Experiments  for  characterizing 
the  spatial  frequency  density  distribution  of  radiant  atmospheric  structure 
(Weiner  spectrum  measurements)  have  been  discussed. 
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As  self-scanning  (CCD,  CID,  etc)  focal  planes  develop  and  extremely 
populous  arrays  become  practical,  real  time  clutter  cancellation  becomes  a 
distinct  possibility.  Dave  Fried  has  suggested  that,  if  detector  spacing  is  at 
least  twice  as  fine  as  the  finest  significant  atmospheric  clutter,  the  clutter 
in  the  track  direction  is  completely  defined  and  its  interference  with  a station 
ary  target  becomes  surmountable.  Weiner  spectrum  information  is  needed 
to  define  the  required  pixel  and  frame  rate  scaling. 
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APPENDIX  B 


TWINKLE 

Observers  at  the  AMOS  facility  began  making  infrared  observations 

of  reentry  vehicles,  artificial  earth  satellites  and  astronomical  objects  in 

1969^^^  with  an  instrument,  designed  at  the  Willow  Run  Laboratory  of  the 

(2) 

University  of  Michigan,  primarily  for  tracking  objects  by  their  thermal 
exitance  in  the  8 to  13  juni  atmospheric  window.  By  1971  the  University  of 
Michigan  Tracker  (UMT)  had  developed  to  the  point  where  satellite  and 
celestial  object  measurements  were  made  routinely  with  it.  Stars  whose 
LWIR  brightness  had  been  measured  and  published  in  the  astronomical 
literature  were  used  as  calibration  standards  for  determining  atmospheric 
extinction  and  the  responsivity  of  our  sensor.  This  author  and  other  ob- 
servers were  puzzled  by  unexpectedly  ragged  LWIR  'calibration  data' 
recorded  even  when  bright  stars  were  accurately  tracked.  As  Figure  B-1 
shows,  large  satellites  on  the  other  hand  occasionally  produced  smooth 
signal  histories  although  tracking  them  was  usually  much  more  precarious. 
In  the  figure,  the  series  of  pulses  is  produced  in  each  case  by  the  image 
sweeping  across  the  detector  field  of  view  stop  approximately  twelve  times 
per  second.  The  field  stop  contains  a picket  fence  reticle  pattern  which 
periodically  occults  the  image  during  each  sweep  producing  a short  burst 
of  1.4  millisecond  pulses  which  are  averaged  to  result  in  a single  pulse, 
approximately  30  milliseconds  wide,  once  ^^er  sweep. 

AMOS  analog  records  of  infrared  star  observations  generally  contain 
random  amplitude  fluctuations  three  to  ten  times  more  prominent  in  the 
8-13  jjjm.  atmospheric  transmission  window  than  in  the  visible.  We  refer  to 
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Figure  B-1  Pegasus  IR  Twinkle 


observations  made  with  the  same  size  (48-inch)  collecting  aperture;  signal 
history  features  which  endure  from  ten  to  several  hundred  milliseconds  and 
detectors  much  larger  than  the  image  diameter  for  both  wavelengths.  Be- 
cause the  IR  intensity  fluctuations  (actually  voltage  fluctuations)  were  much 
more  prominent  than  baseline  noise  and  appear  to  increase  monotonically 
with  atmospheric  path  length  more  rapidly  than  attenuation  scaling  will 
account  for,  they  are  attributed  to  scintillation  produced  by  atmospheric 
refraction  structure  even  though  theory  predicts  little  or  no  difference 
should  be  observed  between  0.  5 and  10  fxm.  Figure  B-1  is  a chart  recording 
of  the  processed  output  of  a UMT  detector  made  during  observation  of  the 
infrared  object  NML  Cyg  (upper  trace)  and  the  Pegasus  U.  S.  satellite, 

SDC  #1384.  The  'star'  is  assumed  to  produce  a X/D  point  image  while 
Pegasus  subtends  over  eight  arcseconds;  strikingly  different  voltage  fluctu- 
ations are  evident.  Observations  similar  to  those  typified  by  Figure  B-1 
led  to  the  notion  that  object  size  information  might  be  embedded  in  endo- 
atmospheric  irradiance  fluctuation  statistics.  Determining  whether  the 
observed  radiation  arises  from  the  object  as  a whole  or  largely  from  some 
small  portion  is  of  interest. 

Actually,  the  observed  fluctuations  may  be  attributed  to  any  of  the 
following  causes:  (a)  system  noise;  (b)  line-of-sight  wander;  or  (c) 
atmospheric  propagation  effects.  In  addition,  the  difference  between  earth 
satellite  and  stellar  variance  may  be  due  to  the  signal  conditioning  equip- 
ment which  includes  an  integrator  to  smooth  out  high  speed  signal  changes; 
the  line-of-sight  moves  100  to  300  times  sidereal  rate  for  near-in  satellites. 

We  have  concluded,  however,  from  observational  evidence  summarized 
here,  that  the  fluctuations  are  probably  due  to  'twinkle'  that  is  to  say: 
distant  refraction  features  in  the  turbulent  atmosphere  producing  interference 
between  light  rays  arriving  at  the  telescope  via  slightly  different  paths. 

With  regard  to  (a)  noise,  note  in  Figure  B-1  that  although  the  star 
(upper)  trace  is  written  with  twice  the  deflection  gain  of  the  Pegasus  record, 
this  in  itself  cannot  account  for  the  smoothness  of  the  latter.  Also,  in  the 
mistrack  region  at  21.8  seconds,  the  background  noise  is  lost  in  the  base- 
line. It  is  therefore  evident  that  peak  signal  fluctuations  are  not  due  to 
instrument  noise  or  radiant  atmosphere  interfering  with  target  flux.  Moreover, 
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it  is  easily  shown  that  satellite  photon  exitance  statistics  cannot  account  for 

such  large  variations.  Figure  B-Z  for  example  relates  noise-in-signal  to 

target  size  and  range  for  surfaces  in  equilibrium  with  insolation  and  earth- 

2 

shine.  For  Pegasus,  with  CA  » 30  meters  observed  at  400  km,  noise-in- 
signal  produces  a negligible  fraction  of  the  observed  fluctuation.  Finally, 
if  system  and  background  noise  are  important  and  relatively  constant  con- 
tributors, normalized  fluctuation  ought  to  be  inversely  proportional  to  target 
irradiance.  This,  however,  does  not  appear  to  be  the  case  as  illustrated 
in  Figure  B-3.  In  the  figure  normalized  standard  deviation,  M,  measured 
during  observations  of  a number  of  infrared  stars,  is  related  to  their  known 
brightness;  there  is  not  enough  correlation  exhibited  to  confidently  tie  the 
two  parameters  together. 

Line -of -sight  wander  or  random  mistracking  convolved  with  detector 

sensitivity  contours,  might  also  account  for  fluctuation  but  this  is  unlikely 

with  UMT.  The  tracker  is  fitted  with  field  lenses  to  provide  a 'Maxwellian 

field-of-view'  for  each  detector.  These  lenses  image  the  telescope  entrance 

pupil  onto  the  semiconductor  surface  such  that,  in  principle,  the  position 

and  intensity  profile  of  the  pupil  image  is  constant  and  independent  of  the 

target  image  within  the  field-of-view  stop.  While  principle  was  imperfectly 

reduced  to  practice'ln  this  case,  ten  arcseconds  of  wander  would  be 
(3) 

required  to  account  for  the  observed  stellar  twinkle.  Ten  arcseconds 
wander  might  be  produced  by  disturbance  torques  due  to  wind  gusts  buffeting 
the  telescope  but  winds  in  excess  of  25  knots  would  be  required  and  observa- 
tions are  rarely  attempted  under  such  conditions. 

Twinkle  strength  should  be  related  to  atmospheric  path  length  and 
Figure  B-4,  which  ties  M%  to  zenith  distance  for  various  stars  observed 
with  UMT,  shows  just  such  relation.  Atmospheric  transmission  dependence 
is  not  severe  enough  (t  = exp  (-0.  1 m),  where  T is  8-13  fjtn  transmission 
and  m is  air  mass)  to  account  for  the  observed  trending.  Unfortunately, 
weather  information  is  unavailable  to  factor  into  these  observations,  note, 
for  example,  the  NML  Cyg  measures  for  28  February  shown  in  the  figure 


3.  Loc  cit  1,  p.  117 
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illustrate  that  weather  changes  may  override  path  length  dependence  on  an 

hour  by  hour  basis.  Generally  February  and  March  are  not  auspicious 

months  for  observations  as  illustrated  in  Figure  B-5.  [Small  circles  in 

the  figure  represent  Zirkind's^^^  manipulation  of  1962-1963  data  obtained 

by  Kuiper  and  associates  at  Haleakala.  The  fact  that  these  values  are  now 

(5) 

believed  to  be  optimistic'  ' does  not  detract  from  the  relevance  of  the 
figure.  ] 

Finally,  a successful  attempt  was  made  to  correlate  twinkle  strength 
with  apparent  angular  dimensions  of  a satellite  tracked  at  high  speed. 

Salyut  #1  (SDC  #05160)  LWIR  data  obtained  in  July  1971  was  used  because 
thij  object  tumbled  end  for  end  about  an  axis  normal  to  the  line-of-sight 
during  the  observations;  it  was  sunlit  and  large  enough  to  produce  resolved 
images  for  a cine  film  record  made  at  the  same  time.  In  Figure  B-6,  the 
encircled  points  represent  rms  deviation  expressed  in  percent  of  the  mean 
signal  value.  Small  photographs  from  the  cine  record,  properly  time 
registered  with  the  IR  scintillation  data,  show  that  target  size  and  not  track 
rate  probably  accounts  for  the  difference  between  the  two  traces  in  Figure 
B-1.  Large  (numbered)  divisions  in  the  scaling  photograph  are  one  milli- 
meter apart  and  the  plate  scale  for  the  60-inch  telescope,  used  to  obtain 
the  cine  recordp  was  8.4  arcseconds/millimeter . Note,  byway  of  example, 
that  Salyut  subtends  almost  eight  arcseconds  at  05:41;  00  U.  T.  The  solid 
curves  in  the  figure,  which  represents  baseline  noise  to  signal  ratio  ex- 
pressed in  percent,  clearly  illustrates  that  signal  fluctuation  is  not  attributable 
to  amplifier  or  detector  G-R  noise.  Data  from  this  and  other  satellite  ob- 
servations were  combined  to  produce  the  dashed  curve  in  Figure  B-7  which 
relates  twinkle  strength  to  object  size. 
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4.  Zirkind,  R.  , Project  AMOS;  An  Infrared  Observatory.  Ap.  Opt.  4,  #9 
p.  1077  (September  1965). 

5.  Miller,  M.  and  Kellen,  P. , Turbulence  Characterization  and  Control, 
RADC-TR-75-185,  p.  5 (July  1975). 
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i Conclusions 


Ground-based  infrared  radiometers  are  size  constrained  and  im- 
practical to  configure  for  imaging  of  resolve  tactical  satellites  but  may 
nevertheless  furnish  gross  size  information  with  intelligence  potential  even 
in  the  absence  of  other  apriori  information  about  the  target.  Size  information 
is  apparently  embedded  in  the  random  variation  about  the  mean  value  of  the 
measured  irradiance  for  very  bright  objects  because  such  fluctuations  are 
produced  by  refractive  atmospheric  eddies  about  the  line-of- sight.  For 
the  8-13  ^m  band  even  a 48 -inch  aperture  is  not  large  enough  to  affect  much 
averaging  and  point  objects  exhibit  appreciable  scintillations  while  objects 
larger  than  two  arcseconds,  which  encompass  several  times  Xj^/r^(ir), 
produce  measurably  less. 

AMT  A and  Twinkle 

The  evidence  presented  so  far  regarding  the  strength  of  LWIR  twinkle 
was  collected  with  the  UMT  which  was  replaced  with  a multicolor  radiometer 
of  higher  sensitivity  (AMTA)  in  1973.  While,  owing  to  press  of  other  business, 
measurements  specifically  designed  to  confirm  UMT  twinkle  results  have  not 
been  conducted  with  the  new  sensor,  there  is  some  evidence  embedded  in 
calibration  exercises.  Figure  B-8  relates  measured  standard  deviation 
about  measured  irradiance  to  star  brightness  for  a number  of  infrared  stars. 
Because  background  flux  exceeds  star  flux  by  four  decades  even  for  the 
brightest  object  observed  here,  s.  d.  , should  not  be  a function  of  brightness, 
yet  we  see  that  it  is.  Stars  above  a threshold  level  established  by  background 
statistics  show  fluctuations  proportional  to  brightness  exactly  as  if  some 
modulating  process-like  line-of-sight  wander  or  atmospheric  optical  effects 
were  operating  on  the  radiation  received. 
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Figure  B-8  NEFD  Related  to  Object  Brightness 
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